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A recently developed scheme [S. Scheel, L. Knòll, and D.- 
G. Welsch, Phys. Rev. A 58, 700 (1998)] for quantizing the 
macroscopic electromagnètic field in linear dispersive and ab- 
sorbing dielectrics satisfying the Kramers-Kronig relations is 
used to derive the quantum local-field correction for the Stan- 
dard virtual-sphere-cavity model. The elèctric and magnètic 
local-field operators are shown to be consistent with QED only 
if the polarization noise is fully taken into account. It is shown 
that the polarization fluctuations in the local field can dra- 
matically change the spontaneous decay rate, compared with 
the familiar result obtained from the classical local-field cor- 
rection. In particular, the spontaneous emission rate strongly 
depends on the radius of the local-field virtual cavity. 

42.50.-p,42.50.Ct,42.50.Lc 



I. INTRODUCTION 

Spontaneous emission by an excited atom is one of the 
most studied examples of a quantum process and may be 
attributed, at least in part, to fluctuations in the electro- 
magnètic vacuum Q. The vacuum field is modified by 
the local environment and this, in turn, leads to a mod- 
ification of the spontaneous emission rate. In this way 
the spontaneous emission rate can be changed by em- 
bedding the radiating atom inside a dielectric host |^| Q , 
or by changing the boundary conditions cither by a cav- 
ity pll·ll^l or a suitable surface [|Ï3|,[Ï4j. Recent experi- 
ments have examined the emission by àtoms embedded 
in dielectric hosts JÏ5|H17[ and have encouraged us to re- 
examine the problem of local-field corrections to the bulk 
modification of the spontaneous decay rate. 

The total decay rate T can be split into two parts, 



r = r J 



m, 



(i) 



in which we associate the transverse decay rate T 1 - and 
the longitudinal decay rate r" with the contributions of 
the transverse and longitudinal fields respectively. The 
dielectric-induced modification of the spontaneous emis- 
sion rate in free space can be ascribed to two effects as- 
sociated with the bulk (macroscopic) field in the medium 
and the other arising from the local (microscopic) field. 
The bulk field correction multiplies the rate by the re- 
fractive index at the transition frequency |||]]. Local- 
field corrections present more of a problem and have a 



form that is strongly model-dependent. For the Clausius- 
Mosotti model, which introduces a virtual cavity sur- 
rounding the atom, a classical treatment of the local-field 
corrections leads to the form MM 



T 



(2) 



for the transverse decay rate of an atom in a bulk di- 
electric of rcfractive index n(u))= y e(ui) — r)(u>) +ík(u>). 
In Eq. (||), r =uJ A fi 2 / (3TTc 3 he a ) is the free-space sponta- 
neous emission rate, where uja and n are, respectively, the 
atòmic transition frequency and the dipole transition ma- 
teix element. The local field correction in Eq. (||) arises 
from writing the local elèctric field in terms of the macro- 
scopic elèctric field and the commonly used induced po- 
larization field. It does not, however, take account of the 
fluctuating component of polarization associated with ab- 
sorption losses. In this paper we investigate the changes 
that arise within the Clausius-Mosotti model when this 
fluctuating component is included. 

Recently, a scheme for quantizing the electromagnètic 
field in an arbitrary linear dielectric medium has been 
proven to be consistent with QED jÏ8|]. It relies on 
the introduction of an appropriately chosen infinite set 
of basic-field operators JÏ9|~pl]] and their connection to 
the electromagnètic field operators via the classical Green 
function. This scheme is a generalization of the approach 
introduced by Huttner and Barnett ^] based on a Hop- 
field model fl23fl of a homogeneous dielectric using Fano 
diagonalization to obtain collective (polariton) exci- 
tations of the electromagnètic field, the polarization and 
the reservoir. As we will show below, the scheme can also 
be applied to local-field quantization that is consistent 
with QED and used for calculating the modified sponta- 
neous decay rate including the local quantum effects of 
the fluctuating medium polarization. 

The paper is organized as follows. After a short review 
of the quantization scheme in Sect. O w e introduce the 
quantum local-field correction in SectTp|. We then apply 
the scheme to the calculation of the spontaneous decay 
rates in Sect. |^ followed by some concluding remarks in 
Sect. [v| Details of the calculation will be given in the 
Appendix. 
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II. QUANTIZATION SCHEME 

We begin with a brief review of the quantization 
scheme used thr oug hout the paper. Further details can 
be found in [[Ï8 20 1. The spectral decomposition of the 
elèctric and magnètic field operators is given by 



E(r) 



(iwE(r,w) + H.c. 



B(r) = / dwB(r,w)+H.c, 



(3) 
(4) 



where É(r, ui) and B(r, lo) satisfy Maxwell's equations 

V-B(r,o;) = 0, (5) 
V- [eoe(r,w)È(r,u;)l =/S(r,w), (6) 



V x E(r,w) = iwB(r,w), 

V x B(r, w) = -i4 e ( r ' w ï( r . w ) + Moj(r, w) 



(7) 
(8) 



[e(r, oj) = eR(r,u) + ie/(r,cj) is the permittivity] . The 
operator noise current density j(r,uj) and the operator 
noise charge density p(r,uj), which had to be introduced 
in order to be consistent with the dissipation-fluctuation 
theorem, are related to the noise polarization P_ N (r,ui) 
as 



j(r,w) =- lW P w (r, w ), 

and satisfy the equation of continuity 
V • j(r,w) = iwp(r,w). 



(9) 
(10) 

(11) 



The operator noise current density j(r,c<j) is obtaincd 
from a bosonic vector field i(r,Lu), 



j(r,u/) = w 



he 



e/(r,w)f(r,w), 



/<(r,w),/t(r>') = dijSir-v'Mu-uj'), 



0. 



(12) 

(13) 
(14) 



The quantization scheme implies that all electromag- 
nètic field operators can be expressed in terms of the 
bàsic fields f(r,w), which may be regarded as being the 
collcctive excitations of the electromagnètic field, the 
medium polarization, and the reservoir. For example, 
the electric-field operator E(r, lo) satisfies the partial dif- 
ferential equation 



V x V x É(r,w) 



such that 



-e(r, w)È(r, u>) = i/x wj(r,w), (15) 



d à swGij (r.s.w^^s^), (16) 



where GV,-(r, s, lü) is the tensor- valued Green function of 
the classical problem. It can then be proven jTJ| that 
this quantization scheme is fully consistent with QED 
for arbitrary linear dielectrics. 



III. QUANTUM LOCAL-FIELD CORRECTION 

If we think of an atom located at some space point 
inside the dielectric, then the macroscopic field of Sect.|| 
will not, in fact, be the field felt by the atom. From clas- 
sical electrodynamics we know that we should introduce 
what is called the local field at the location of the atom 
p5j. There are essentially two ways of introducing the 
local field. First, one could cut out a real cavity (most 
commonly a sphere) around the atom and calculate, in 
our scheme, the elèctric field inside the cavity according 
to Eq. (|Ï6|). This would lead us to introduce the elèctric 
field operator E loc (r,o;) by the relation 



^ oc (r, W )= W 



d 3 sc^f(r, S)W )j.(s,u;), 



(17) 



where GJ" h (r,s,w) is the Green function of the classical 
problem of an inhomogeneous medium that consists of 
the real cavity surrounded by the dielectric in which the 
atom is embedded. 

To avoid the solution of the inhomogeneous problem, 
commonly a simpler virtual-c&vity model of Clausius- 
Mosotti-type is used. In this model the local field is (ap- 
proximately) related to the macroscopic field, which in 
our scheme can be obtained according Eq. ( |Ï6| ) , with the 
Green function for the bulk-medium problem. In classi- 
cal òptics the elèctric field is usually corrected by adding 
a term proportional to the polarization, 



where 



E loc (r,c) =E(r,^) + — P(r,u), 
3e 



P(r,^) = 6 [e(r,w)-l]E(r,w), 



from which it follows that 



E loc (r,^) = i[6(r, W ) + 2]E(r^). 



(18) 



(19) 



(20) 



This classically corrected local field is just the field used 
for the derivation of the rate formula (||) . 

In classical òptics Eq. ( [ÏTj| ) corresponds to the zero- 
temperature limit. It does not, however, hold in opera- 
tor form in quantum òptics. Rcgarding Eq. (19) as an 
operator equation would strongly contradict the QED 
equal-time commutation relations in general, so that the 
operator counterpart of Eq. ( p0| ) will not be a true elèctric 
field. In order to obtain a canonical operator equation, 
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we have to complement Eq. ( |Ï9| ) by a noise polarization 
P N (r,u): 



P(r, w) = e [e(r, u) - l]E(r, a,) + P w (r, w 



(21) 



Using the quantization scheme in Sec. ||, from Eq. (g) 
the noise polarization P (r, u;) is given by 



P N (v, lj) = -- j(r, c) = iJ ^ ej (r, W ) f (r, «). (22) 

- V 7T 



Combining Eqs. (|18|) and (|21|), a local electric-field oper- 
ator 

È loc (r, W ) = § [e(r, w ) + 2] È(r, w) + ^- P K (r, w) (23) 

àeo 



can be introduced, withP Ar (r,w) being given in Eq. (22). 

In order to prové the consistency of the quantized local 
field with QED, we compute the (equal-time) commuta- 
tion relation between the fundamental local fields E loc (r) 
and B loc (r). For this purpose we note that the local 
magnètic- field operator is given by, on using Eqs. ([?]), 



B loc (r,cj) = V x V— È loc (r,w), 



(24) 



where the symbol V stands for the principal part. Re- 
calling Eqs. (||) and (0), the local-field operators in real 
space are 



/>oo 

È loc (r) = / du}É loc (r,uj) +H, 
Jo 



and 



B loc (r) 



dwB loc (r,cj) +H.c. 



(25) 



(26) 



Expressing the local elèctric and magnètic fields in terms 
of the bàsic fields f(r, u), from the calculation given in 
App. |À| it is found that 



Èl oc (r),Èt c (r') 

ih 



0) 



eiW ^(r-r'){l + i[e(r,0)-l]}, (27) 



and it is easily seen that 

"£i oc (r),£ fc loc (r')l - [^(r),^-^') 



= 0. 



(28) 



The result reveals that the local-electric field operator 
( p3| ) and the associated magnètic field operator ( p4| ) can 
be regarded as being consistent with quantum theory, 
provided that the (real) static permittivity es(r)=e(r, 0) 
satisfies the condition 



es(r) 
10 



< 1. 



(29) 



Equivalently, the static refractive index ng(r) = y'es(r) 
must be small compared with vT0«3.16. 

It should be noted that a term proportional to the 
5 function 5(w) can be added to the right-hand side of 
Eq. (|24|) in order to recover Ampere's law when the equa- 
tion is multiplied by lj. Obviously, this ambiguity reflects 
the fact that the static magnètic field cannot be inferred 
from the static elèctric field. From a simple calculation it 
can be shown that such a term does not change the com- 
mutation relation (|2^). Since it is only relevant at zero 
frequency, it does not play any role in the calculation of 
the decay at transition frequency u>a- 

In order to take into account a possible deviation of 
the symmetry of the material from cúbic symmetry, a 
structure constant s can be included in Eq. (|ï^) such 
that l2l 



E loc (r,^) =E(r,oj) + - 



1 

3 +S 



P(r,w). (30) 



Regarding this equation as an operator equation with 
P(r,oj) from Eq. ( pï|) and following the line in Ap- 
pendix it can be seen that Eq. (|2^) changes to 



£] oc (r),i?,! oc (r') 



= -^e ife ^(r-r')|l+y[e(r,0)-l]|, (31) 

where the parameter a in related to s by 

a = 1 + 3s. (32) 

Thus, consistency with quantum theory is achieved, if 
the condition 



es(r) < 9a 2 + 1 



(33) 



is fulfilled. 



IV. SPONTANEOUS DECAY RATE 

The spontaneous decay rate of a (two-level) atom at 
transition frequency uja placed at point ya is given by 

r = — J duj fj,i{0\E} oc (r, u)E_] oc] {v A , u; A )\0)fJ,j (34) 

(r — > ta)- In what follows we consider a homogeneous 
bulk material, i.e., e(r, u>) = e(u>), and assume that the 
inequality ( p9| ) is fulfilled. Using Eq. (p3|), the vacuum 
expectation value of the local electric-field operators in 
the limit r — ► can be written as the sum 

(011^(^)1^(^)10) 



e(w) + 2e*(w') + 2 



|4(r,c)l/(r A , W ')|0) 



+ 



$4 



\È?(t,v)ÈF(ta,«/)\0) 



3 



eH + 2 

9e 
e*(iü') + : 

9e 



(0\P?(r, U )Ètir A ,u>')\0), 



(35) 



where É(r, uj) and P_ N (r,oj) are given by Eqs. (|ï|) and 
( f22| ) respectively. The Green function G^(r, r^, uj a ) for 
the bulk material in the limit r — * has the form 
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Gij (R, oj A ) = Gfj (R, u A ) + G% (R, u> A ) 
(R = r — ta), where 



(36) 



GtAK,Lü A ) = — 



1 í RiRj 



4tt 2iï 3 



2i? 



2zo;,4 
3c 



[t)(uj a ) + ík(oj a )]Sí 



O(R) 



(37) 



and 

Gl(K,w A ) = - 
iir 



AiïU) 2 A e{w a) 



3 RiRj 



1 

R? 



(38) 



are the transverse and longitudinal parts respectively. 
We see that the real part of the transverse Green function 
as well as the longitudinal part itself diverge as R— *0, 
reflccting the fact that a macroscopic approach is vàlid 
only to some appropriately fixed scale R which exceeds 
the average distance of two àtoms in the dielectric. Fol- 
lowing [Q, we average the divergent terms in the decay 
rate over a sphere of radius R. We will see later that one 
can estimate an upper bound for R from the decay rate 
itself. 

The first term on the right-hand side in Eq. (|35| ) gives 
the contribution to the decay rate with the classically 
corrected local ficld Mm, 



r r . 



. r 



e{w A )+2 


2 


3 





2\e{Lü A )\ 2 \cü A R 



(39) 



with the transverse rate T^j being given in Eq. (j^). The 
second term in Eq. ( (35|) is purely a contribution of the 
noise polarization held and is given by 

1 



^(0\P?(t,u)PP{t A ,ü/)\0) 



9vre 



SijSÇR) 6{u - ui'). 



(40) 



The cross-terms mixing the macroscopic elèctric field and 
the noise polarization field give rise to the contribution 

e(ui) + 2 



+ 



(jío (0|4(r, W )£ft(rA,</)|0) 



9e 

2w 2 h , , 
— — e/(w)Rc 



Pf(r,u)EUr A ,u')\0) 



3ttc 2 eq 



e(Lü) + 2 



5{üj-üj'). (41) 



Hence, the total decay rate reads 
2 



r c i + 



9& 

4 «4 



ei(u) A )iMHj Sij6(R) 

e(iü A ) + 2 



3Heo c 2 



ej(u A )fJ,HJíjBe 



(42) 



Equation (|42| ) is remarkable in the sense that inclusion 
of the polarization noise in the local field gives rise to 
a term that only results from that noise and leads to 
a dependence of the decay rate on the real part of the 
Green function. We average the 6 tensor 

^•5(R)=^(R) + 4(R) (43) 

and the Green tensor ( p6| ) over the sphere and obtain 

1 



^(RH2^(R) = _%, 

1 W A k((V A ) 



ReG+(R,au) 



ImG|.(R,a; A ) 



inR 6irc 

c 2 e R (üj A ) , 

47r W 2| e ( WA )|2 j R3 
ÜJ A T)(W A ) 

6irc 131 
c 2 €i(uja) - 
±™ A \e(u> A )\ 2 R 3 ij ' 



(44) 

(45) 
(46) 
(47) 
(48) 



and Eq. fl42|) can be given in the form of Eq. (@), where 
T- 1 - and rllread as 



e(üJ A ) + 2 



-ei(u> A ) [e R {u A ) + 2} 
c?ei(u) A ) \ 



2e 2 {üj A ) 



2 t \ 



3ujaR 



+ 



3uj\R 3 



(49) 



and 

r ll = 



^ 2 ei{ujA) 



2nhe R 3 \e(uA)\ 2 
2 



e{u> A ) + 2 



<o>a)\ 1 



9 



e R {aj A ) [e R (uj A )+2] - - [ei(u) A )] 



(50) 



The modifications near a medium resonance are clear. 
Note that owing to the quantum local-field correction the 
unspecified parameter R also enters into the transverse 
decay rate. In order to compare our canonical result with 
that obtained using the classically corrected local field, 
we use, for comparison, the same Lorentz model Jor the 
permittivity of a single-resonance medium as in 



4 



e(o;) 



(0.46 lü t ) 2 

\ — LO 2 — ÍJLÜ 



(51) 



where lot is the resonance frequency of the medium. Fig- 
ures show the transverse decay rate T 1 - with and 
without quantum local-field corrections as a function of 
the atòmic transition frequency loa for different vàlues 
of the damping parameter of the medium, 7, and the 
parameter 



(52) 



First of all, for small r, i.e., large virtual-cavity radius 
R one observes little reduction of spontaneous decay for 
freqüències loa just above the resonance frequency lot- 
Its possible applications in semiconductor physics and 
solid-state physics has already been discussed J27J . 

The greatest difference between the quantum-mecha- 
nically and classically corrected transverse decay rates 
T 1 - and Igj, respectively, arises near the medium reso- 
nance when 7 is small. Both the imaginary part of the 
permittivity and the real part can take very large vàl- 
ues for lo a ~ wr and in consequence T 1 - can drastically 
change compared with T^. Obviously, in the resonance 
regime the noise polarization essentially contributes to 
the local-field and therefore strongly influences T^. For 
small vàlues of 7 both qualitative and quantitative diffcr- 
ences between the rates T ± and T^j are observed (Fig. [!]). 
With increasing value of 7 the two rates become less dif- 
ferent from each other, the changes being quantitative 
rather than qualitatitive (compare Fig. [ï] with Fig. [|) . 

Whereas is independent of the radius R of the vir- 
tual cavity, T -1 - sensitively depends on R. To fix the value 
of R that is undetermined in the Clausius-Mosotti model, 
experimental data could be used in principio (for recent 
experiments on spontaneous emission, see, e.g., JÏ5|-p7t). 
It is worth noting that our theory gives an upper bound 



Rn 



i.e., a lower bound r m ; n for the parameter r, be- 



cause of the fact that r cannot be negative. As already 
mentioned, the limit r — > ya in Eq. (p3) cannot be per- 
formed and averaging over a sphere of radius R can give 
negative vàlues, if the radius is not small enough, because 
the vacuum expectation value of the Hcrmitian operator 
É loc (r, w)È loc '''(r J 4, lo) is not necessarily positive. Figure^j 
presents r min as a function of the damping parameter 7. 
The curve was obtained numerically by requiring that T 1 - 
must not be negative over the whole frequency spectrum. 

Figure shows that for chosen (small) 7 and r < 
r m j n negative vàlues of r may appear when the atòmic 
transition frequency loa approaches the medium reso- 
nance frequency lot and is in an interval that corre- 
sponds to the polariton bandgap between lot and lol — 
[lo^+^A&lot) 2 ] 1 ^ 2 in the Hopfield model of a dielectric in 
the absence of absorption |[23| . Obviously, in this regime 
of spontaneous decay, in which the decay rate is expected 
to respond to a change in the local-field noise very sen- 
sitively, the Clausius-Mosotti model has to be replaced 
with a refined model, at least in quantum theory. 



From the Standard derivation of the (classical) Clau- 
sius-Mosotti local-field (see, e.g., the radius R of the 
virtual cavity should be larger than the average distance 
of two neighboring àtoms but sufficicntly smallcr than 
the optical wavelength. In terms of the parameter r, the 
latter requirement means that r^$>l. Provided that the 
damping parameter 7 is not too small, this is in agree- 
ment with the condition that the parameter r should not 
be smaller than r m ; n given in Fig. ||. 



V. CONCLUSIONS 

Within the frame of the Clausius-Mosotti model we 
have studied the influence of the quantum local-field cor- 
rection arising from the noise polarization on the spon- 
taneous decay rate of an excited atom embedded in an 
absorbing medium. We have shown that inclusion in the 
local field of the noise polarization ensures that the lo- 
cal field fulfills the fundamental equal-time commutation 
relations of QED, provided that the static refractive Ín- 
dex of the medium does not exceed unity substantially. 
The calculated rates demonstrate that the contribution 
of the noise polarization to the local field is extremely 
important and cannot be ignored. In particular, at the 
resonance freqüències of the medium the transverse decay 
rate can drastically change compared with the classically 
corrected rate where the fluctuating component of the 
polarization is omitted. 

The decay rate crucially depends on the choice of the 
radius of the virtual cavity. It is worth noting that from 
the dependence on the radius of the transverse rate - an 
effect that only results from the noise polarization - a sec- 
ond condition of validity can be imposed on the local-field 
model. In order to obtain for any transition frequency a 
positive transverse decay rate, the cavity radius must not 
exceed some upper bound. 

The Clausius-Mosotti virtual-cavity model is com- 
monly based on the assumption that the near field that 
arises from the àtoms inside the cavity averages to zero 
in general. In quantum òptics this assumption may fail, 
because of the vacuum noise associated with this field, 
which may be an explanation for the restrictions found. 
In order to extend the range of validity of the theory, a 
more refined concept seems to be necessary that consid- 
ers the near field more carefully. 



5 



APPENDIX A: COMMUTATION RELATIONS OF THE LOCAL FIELD OPERATORS 

From Eq. (|2^) together with Eq. (^), the local elèctric field operator reads in Fourier space 

e(r,w) + 2 



E loc (r,w) 



•-! : - (r ' w) + 3^V í £/(r ' w) ^ r ' w) 



Combining Eqs. (|24|) and (Al), we obtain for the local magnètic field in Fourier space 
B loc (r', W ')-Vx 



1 . (s e(r / ,a/) + 2 1 1 /fien , , . 
iu/ 3 3en w y i 



(Al) 



(A2) 



Recalling Eqs. ( p5| ) and (pq), from Eqs. (Al) and ( |A2| ) together with Eq. ( |Ï6| ) the local elèctric and magnètic fields 
are given by 



e(r, w) + 2 íüü 



3e 



— ^/ e/(s, w) (r, s, w)/ i (s, w) + H. 



f v-, 

3e \ n J Q 
fien 



B fc loc (r') = e Wm df {\I'-^V I du' / dV 



OO 



1 /^o^ 
3e V 7r j 



o 



£(r, ' Q U;,) + 2 ^ >/e/(s', w') G m „(r', s', u/)/ n (s', (/) + H.c. 
3e c z 



V^y)/m(r'w')+H.c. 



Thus, the (equal-time) commutator between the local elèctric and magnètic fields can be given by 
£>° c (r),B]fV)" 

~ílü 3 e(r,w) + 2 e*(r',w) + 2 



fie, 



K e klm dUVj du j 



3e 



3e 



e / (s,w)G íj (r,s,w)G* ni (r',s,w) - c.c 



1 






V 


3e 




1 


V 


3e 




2i 






V 


94 





du 



e*(r',u;) + 2 íüj 
3e c 2 



ei(r,Lü)G* mi (r',r,uj) - c.c. 



e(r,w) + 2 iu , , 

e/(r , u)G im {r, r , w) - c.c. 

3e c z 



-- [e/(r,^ ro í(r-r')] 



The remaining spatial integral in Eq. (A5) can be calculated using the symmetry relation 

G tj (r,r',uj) = G^r' ,r,u), 

the crossing relation 

G 4j ·(p,r' 1 w)=GÍ.(r 1 r',- W ), 

and the integral relation p(| 

/ d 3 se/(s,w)G K (s,r,w)Gy(s,r',w) = — [G^r', r, w) - G?-(r, r', u)] 
Straightforward calculation yields 



7re 



/£>(r,r')+J^(r,r') 



(A3) 



(A4) 



(A5) 

(A6) 
(A7) 

(A8) 
(A9) 



where 



l£>(r,r') = V I ^G iro (r,r',u;) {l + ± [e(r,u;)-l] + § [e(r»-l] + i [e(r, W )-l] [e(r',u) - 1]} , (A10) 



e/(r,w) 



<5im<5(r - r'). 



(All) 



Closing the integration contour in the upper complex frequency half-plane and following the line in lla ], we derive 
that 



lM(r,r')=inS im 5(r-r'). 



Recalling the Kramers-Kronig relations, the lj integral in Eq. (All) is easily performed to obtain 

/j£ (r, r') = tt [e R (r, 0) - 1] S m S(r - r') = tt [e(r, 0) - 1] 5 im 5(r - r'). 
Combining Eqs. (A.9), (A12), and ( A13| ) then yields the commutation relation 



(A12) 



(A13) 
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FIG. 1. The transverse decay rate T x is shown as a func- 
tion of the transition frequency loa for 7 = 0.01 and r = 10 
(dashed curve), r = 20 (broken curve) and r = 30 (dotted 
curve). For comparison, the rate without quantum local- field 
correction j| is shown (sòlid curve). Since for r = 10 (broken 
curve) F x becomes negative, this case must be excluded from 
consideration (cf. Fig. 
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FIG. 2. The transverse decay rate T x is shown as a func- 
tion of the transition frequency loa for 7 = 0.05 and r = 10 
(dashed curve), r — 20 (broken curve), and r = 30 (dotted 
curve). For comparison, the rate without quantum local- field 
correction M is shown (sòlid curve). 
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FIG. 3. The transverse decay rate r is shown as a func- 
tion of the transition frequency uja for 7 = 0.1 and r = 10 
(dashed curve), r — 20 (broken curve), and r — 30 (dotted 
curve). For comparison, the rate without quantum local- field 
correction [H] is shown (sòlid curve). 
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FIG. 4. The lower bound r m i n of the parameter r, 
Eq. (^), is shown as a function of the damping parameter 
7. The region below the curve is the part where the trans- 
verse decay rate F x may take negative vàlues and is therefore 
forbidden. 







